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Abstract 

Within the diffusion Monte Carlo approach, we have determined the structure of isotopically 
pure and mixed hehum droplets doped with one magnesium atom. For pure ^He clusters, our 
results confirm those of M. Mella et al. [J. Chem. Phys. 123, 054328 (2005)] that the impurity 
experiences a transition from a surface to a bulk location as the number of helium atoms in the 
droplet increases. Contrarily, for pure ^He clusters Mg resides in the bulk of the droplet due to 
the smaller surface tension of this isotope. Results for mixed droplets are presented. We have also 
obtained the absorption spectrum of Mg around the 3s3p^Pi ^ 3s^ ^5*0 transition. 

PACS numbers: 36.40.-c, 33.20.Kf, 67.60.gj 
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I. INTRODUCTION 



The solubility of alkaline earth atoms attached to helium droplets has been found to 
critically depend on the dopant and on the helium isotope droplets are made of. This is at 
clear variance with what happens for alkali atoms, that are all found to reside in a dimple at 
the surface of the droplet irrespective of their isotopic composition,-"^ or to impurities that 
experience a large attractive interaction with helium as, e.g., inert gas atoms that are found 
to reside in the bulk of the droplet. This is a particular yet prominent aspect of a much 
broader subject -the physics and chemistry of pristine and doped helium nanodroplets- that 
has been reviewed in a series of articles, see e.g. Refs. |3|-[8| and references therein. 

Ancilotto and coworkers^ have provided a solvation criterion for impurities in liquid he- 
lium in terms of the dimensionless parameter 

where p and a are the density and the surface tension of liquid He, respectively, and e and Ve 
are the well depth and the equilibrium distance of the He-impurity interaction, respectively. 
This parameter measures the balance between the energy of the impurity and the surface 
energy of the liquid. If A < 1.9, the impurity sits on the free surface of the fluid and no 
solvation occurs. For Mg one has A = 2.6 for ^He and A = 4.6 for '^He. Since this criterion 
does take into account neither the complexities of the system nor the fine details of the 
He-impurity pair potential and it has been established for bulk liquid, it must be taken with 
care when the value of A is fairly close to 1.9, in which case only a detailed calculation may 
unveil the solvation properties of a given impurity in He drops. 

It has been found experimental and theoretically that Ca atoms solvate in ^He but not 
in ^He droplets, and that depending on their isotopic composition, they may reside in the 
^He-^He interface that develops in mixed droplets at low temperatures.—"— The lighter 
alkaline earth Mg presents a borderline behavior. Within Density Functional Theory (DFT), 
it has been found to solvate in ^He droplets and to be very delocalized in ^He ones.— 
Previous diffusion Monte Carlo calculations (DMC)^'^ have yielded the result that Mg is fully 
solvated for a critical number of ^He atoms of about N4 = 30. A similar transition from 
surface to bulk location of Mg was also found within DFT.— Experimentally, full solvation 
of Mg in ^He droplets has been inferred from the analysis of Laser Induced Fluorescence 
(LIF),— comparing it with LIF experiments on the absorption spectrum of Mg atoms in 



liquid ^He.— More recent Resonant Two-Photon-Ionization (R2PI) experiments^^ have 
also yielded a bulk location for this impurity. 

Since an electron-impact ionization experiments^ carried out in Mg doped ''He drops with 
about 10^ atoms seems to indicate that magnesium reside at or near the droplet surface, 
and some unpublished DMC calculations seem also to point toward a rather surface location 
of Mg in ^He droplets,—"— we have undertaken an independent Monte Carlo analysis of the 
structure of these systems aiming at settling this issue, complementing the previous analyses 
with a DMC study of ^He and mixed ^He-'^He droplets doped with Mg that has not been 
previously carried out. The interest in addressing doped mixed helium clusters has been 
already stressed.-ii^ Very recently, we have studied the structure and absorption spectrum 
of Mg in liquid helium mixtures as a function of pressure and isotopic composition.- 

This paper is organized as follows. In Sec. [TTl we give some details about the DMC 
calculations. In Sec. IIIII we present our results, and in Sec. IIVI we present a brief summary. 



II. METHODS 



Quantum Monte Carlo methods aim at solving the Schrodinger equation of a many-body 
system with the only knowledge of the inter-particle interaction. In our calculations we 
have employed the He-He Aziz potential,— and the X^S Mg-He interaction as obtained by 
Hinde.— For the description of the absorption spectrum described in Sec IIII C[ the S and H 



Mg-He pair potentials of Ref. [ij have been used. To facilitate the use of these potentials, 



we have fitted them (by simulated annealing) to an analytical expression of the kind 



VM,-H.{r) = Ae—^- - Fir) + ^ + (2) 



F(r) 



1 r>D 

The parameters are given in Table [B The resulting fits are shown in Fig. [H together with the 
calculated values (indicated with symbols in the figure). It can be seen that the agreement 
is excellent. For the sake of comparison, we have also included the He-He interaction used 
in this work. 

Our DMC calculations are based on a variational or importance sampling wave function. 
In that case, one does not solve the Schrodinger equation for the true wave function \E'(7^, t), 



but for the auxiliary function f(7l,t) = The trial wave function \&t(7?.) 

guides the random walk and minimizes the variance. We have used a rather simple form, 
containing the basic required properties, the same as employed in the past to describe a Ca 
impurity. It is a generalization of the trial function adopted in previous studies on pristine 
mixed helium clusters,— and is written as a product of seven terms 

^y(7^) = ^33^34^4Mg^3Mg^t^i ' (3) 

where {7^} represents the set of 3(Ai'3 + A^4 + 1) coordinates of the atoms forming the cluster, 
^3(^4) being the number of ^He (^He) atoms. The first five terms are Jastrow factors "^mn 
for each pair (M, N) of different atoms, for which we have chosen the generic form 



- OiMNrij , (4) 



where indices i,j run over the corresponding type of atom, and includes a short-range repul- 
sion term with parameters Bmn and i'mn associated with it, and a long-range confining term 
with corresponding parameter omn- The form of the short-range repulsion term was intro- 
duced long ago by McMillan^^ to describe the homogeneous liquid ^He using a 12-6 Lennard- 
Jones interaction. In that case, the values of the two parameters u and b are fixed by the 
short-range behavior of a pair of atoms, and one gets z/ = 5, and b = [16/ie/(25/;.^)]^/^''cr^/^, 
where /i is the reduced mass for each pair of atoms, and e and a are the Lennard- Jones 
energy and distance parameters, respectively. 

Notice that each "^mn function is explicitly symmetric under the exchange of identical 
particles. The antisymmetry required for ^He fermions is incorporated in the Slater deter- 
minants and D^, related to the spin-up and -down fermions, respectively. Two aspects 
should be considered in these determinants, namely the form and the filling of the single 
particle orbitals. As in previous works^^^^^ the determinants have been built up with homo- 
geneous monomials of the fermion Cartesian coordinates as x^"" y^^ z^'' , where the subindex i 
refers to the particle, and the integer number n = n^^ + Uy + defines a "shell". This kind 
of Slater determinant has also been employed for other many-body systems, (see e.g. Ref. 
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29|). It turns out that they are of Vandermonde type and, provided that the shells are 
filled in an increasing order of ra, they can be expressed in terms of products of the relative 
coordinates, being thus translationally invar iant.^--^ Besides, we have always assumed a 
filling scheme in which the total spin is minimum, either or 1/2, for A'"3 even or odd, re- 



spectively. For some pure ^He clusters, the shell filling has been given in Ref. |26|, indicating 
the total value of Sz- 

For the DMC algorithm we have used the short-time Green function approximatiort^i^ 
with an O(r^) form.^ For those systems involving nodal surfaces, i.e. for fermions, the fixed 
node approximation has been employed.— i^i^ The energy accuracy depends on the quality of 
the nodal surfaces of the trial function, which arise from the Slater determinants. Feynman- 
Cohen back-fiow correlations^ have been incorporated into the scheme by substituting 

Vi^ii = ri + Y, - r^) (5) 

in the orbitals of the Slater determinants.— For the backfiow function ri{r) we choose the 
medium-range form used in Ref. |37|, namely ri{r) = A/r^, with the same value of A = 5 A^. 



These type of correlations give rise to nodal surfaces that provide very accurate results 
for different fermionic systems in Quantum Monte Carlo calculations.—'^"— Pure estimator 
results discussed later on are free of the trial wave function bias in the density distributions. 
The accuracy of the wave function only affects the rate of the convergence of the calculation. 

As in our previous works on pure and doped mixed helium clusters,— >2^i^ we have fixed 
the parameters Bmn by taking e as the minimum of the interaction, and a as the distance at 
which the interaction is zero. Moreover, we have slightly modified the value of the exponent 
u with respect to the McMillan value. All in all, for all cluster sizes we have used the 
following values: umn = 5.2, 644 = 2.95 A, 633 = 2.85 A, 634 = 2.90 A, 64Mg = 5.831 A, 
and 63Mg = 5.687 A. Thus, the trial or importance sampling wave function contains only 
five free parameters, namely ^44, 033, ^34, a^ug, and a3Mg, which have been determined 
by minimizing the expectation value of the Hamiltonian. It turns out that for isotopically 

pure boson (A^3 = 0) or fermion (A^4 = 0) clusters, these parameters are well represented 

1/3 

in terms of the number of atoms in the drop, namely 044 = 0.117 — 0.029A^4 , a4Mg = 
0.268 - O.OllNl^^, 033 = 0.020 + O-OOGN^^^, and a3Mg = 0.191 + OmN^^^ - 0.152N3^\ 

Our DMC calculation is based on an importance sampling function, and therefore the 
walkers are generated according to the auxiliary function f{TZ,t) = \1't'(7^)\E'(7^, t). The 
natural output corresponds thus to the so called mixed estimator, in which the expectation 
value of a given operator is straightforwardly calculated with this probability distribution 
function. If such an operator commutes with the Hamiltonian, the mixed estimator is equal 
to the exact expectation value of the observable in the asymptotic limit and within the fixed- 



node error. For other observables, as for example radial operators, the mixed estimator is in 
general biased by the trial function used for importance sampling. For the reasons mentioned 
in the Introduction, we have chosen to obtain unbiased estimates of radial distances for some 
^He droplets. In view of the few unbiased available calculations, ours may contribute to the 
general discussion of the subject. 

Within Quantum Monte Carlo simulations, several schemes have been proposed in the 
literature to obtain unbiased — also called pure — estimators of operators that do not com- 
mute with the Hamiltonian. In this work we compute pure expectation values by using a 
forward walking method. The idea is to recover the exact value by including the factor 
^(7?., t)/\l/r(7?.) in the expectation value. This quotient can be obtained starting from the 
asymptotic offspring of the walker.~ Therefore a weight proportional to the number of the 
future descendant of the walker 

V^(7^) oc r^(7^,^ ->cx)) (6) 
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needs to be included in the calculation. We have employed the algorithm devised in Ref. 
to calculate this weight. It makes use of an auxiliary variable associated with each walker 
that evolves with it, i.e. it is replicated as many times as the walker and propagates the 
local values. With the proper boundary conditions, the final average provides the pure 
expectation value. The basis of the method and details on the algorithm can be found in 



Refs. 
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The pure estimation depends on the size of the block, A^, that needs to be 
large enough to fulfill the forward walking condition Eq. (j6]). A study of the convergence as 
a function of is required to fix the block length from which the pure estimator provides 
the same value within the statistical error for a given radial operator. 



III. RESULTS 



A. Structure of Mg@^HeAr4-|-'^Hejv3 clusters 

For illustrative purposes, the mixed estimator total helium particle densities of the "^He4o, 
^He2o+^He2o and ^He4o droplets doped with Mg are plotted in Fig. O These densities have 
been obtained as follows.^ After a simulation running for a long thermalization time, we 
have stored a large number of walkers (typically 10^) and for each of them the origin is 
taken at the center of mass of the helium atoms. Next, a rotation is carried out so that the 
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Mg atom lies on the z-axis. In this coordinate system, a projection onto the y = plane is 
performed in order to compute the density. 

It can be qualitatively seen that, while for ^He4o Mg is fully immersed, for '^He4o it is not. 
We will see in Sec IIII CI how the distinct helium environment around Mg is reflected in the 
absorption spectrum for these three complexes. 

Figures [3] and H] represent the mixed estimator densities corresponding to the doped 
droplets ^He8+'^He2o and ^He2o+'^He2o, respectively. It is interesting to notice that the ^He 
component may get in touch with the impurity. This is even so for the A^s = A^4 = 20 
system, the reason being that the balance between the weak He-He and Mg-He interactions 
favors that the Mg atom is not fully coated by ^He for such small droplets and the chosen 
composition. The situation changes when A^4 increases.—"^ 

The mixed estimators for the root mean square (rms) radius of the helium cluster and 
the Mg impurity are plotted in Fig. [51 We have chosen to refer them to the center-of-mass 
of the He droplet instead of the center-of-mass of the He-I-Mg complex, so that the data 
clearly show where the Mg sits with respect to the He moiety. 

In ^He droplets the Mg rms radius increases with cluster size for small (A^4 ^ 25) droplets. 
One may infer from this that, for small A^4 values, the impurity sits in the outer region of the 
droplet. At A^4 ^ 25 the impurity begins to "sink" into the cluster as its rms radius starts 
to decrease with increasing cluster size. In ^He droplets no such trend is seen for the Mg 
rms radius, which contrarily has a tendency to decrease as the cluster size increases, apart 
from the structure around N3 = 14 for whose origin, likely related to shell effects, we have 
been unable to find a convincing explanation. These results are consistent with the finding 
that Mg is in the bulk of the droplet for any value,— as Ancilotto's criterion predicts. 

The transition in the Mg location Fig. [5] hints at can be clearly seen in Fig. [6l where the 
mixed estimator radial probability distribution of Mg in ^He clusters is shown for several 
sizes from A^4 = 8 to 30. For small cluster sizes (A^4 < 25), the Mg atom is found at the 
surface of the droplet. The typical distance between the impurity and the center-of-mass of 
the droplet is about 3-4 A. At A^4 = 25 Mg experiences a transition from a surface to a bulk 
state in which the probability of finding the impurity inside the droplet becomes significant. 
For cluster sizes between A^4 = 25 and 29, the impurity is highly delocalized, and for A^4 > 30 
the Mg atom resides inside the droplet. The typical distance between the impurity and the 
center-of-mass of the droplet decreases as the impurity becomes fully solvated. Our results 
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thus confirm those of Mella et al— 

Finally, we have employed the pure estimator to determine the radial structure for the 
isotopically pure A^4 = 8 and 30 droplets. We have found that the differences between pure 
and mixed estimators of the radial densities are very small and not worth to be plotted. As 
an illustration, we have shown in Fig. [7] the mixed and pure estimations for the rms radius 
as a function of the block size Al. The stability of the results is apparent, as well as the 
quality of our mixed estimator for radial distances. Thus, the conclusions drawn using the 
mixed estimator are robust and remain unchanged and accordingly, apart from Fig. [71 the 
results discussed in this work have all been obtained using the mixed estimator. 



B. Energetics of MgO^HeAr^+^HeTVa clusters 

The calculated DMC ground state energies E(Mg@^}ieiy^ +^HeAr3) are given in Tables HTl 
IIIH and [IV] for several (N^jN^) combinations. In these Tables we have shown the statistical 
error in the last figure obtained in the usual way, as the standard deviation of the mean 
value calculated by using the blocking method. The solvation energy of the dopant in the 
droplet, defined as 

/XMg = E{Mg@^ReN, He^3) - ECRbn, He^3) , (7) 

is also given in the Tables. No solvation energy is given for MgtQ^HeATg droplets because pure 
^He droplets are unbound for such small sizes^^'^'^ and hence the solvation energy is just 
the ground state energy of the MgQ^HeAfg complex. At first glance, it is surprising that one 
single Mg atom is able to bound any number of ^He atoms, since the Mg-He interaction is 
weaker than the He-He one, see Fig. [1] The smaller zero-point energy of Mg together with 
the fact that the equilibrium distance is larger for Mg-He than for He-He causes the extra 
binding. Thus, the complexes MgtQ^HeAr^, Mg@^HeAr3 and Mg@^He7V4-|-^HeAr3 are bound for 
any and A^4 values. 

A comparison between our DMC results for MgtQ^HeAr^ clusters and other DMC results 



from Refs. 



14 and 



221 is presented in Table IVl Mella et al}^ have used a Mg-He interaction 
determined at the CCSDT level, and the He-He interaction of Ref. 0. The differences 
beyond statistical errors found between our results and theirs could be mostly attributed 
to minor differences in the interaction potentials. The agreement with Elhyiani^S -who has 
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used the same interactions as in this work- is very satisfactory. 

The ground state and Mg solvation energies are presented in Figs. [8] and [9] respectively, 
as a function of the total number of He atoms. For the sizes considered here, it can be 
seen that the Mg solvation energy does not much change when ^He atoms are added to the 
Mg@4He^4 droplet. 

Both ground state and solvation energies are smooth functions of the total number of 
helium atoms. It is worthwhile noting that for doped '^He droplets, these energies display 
conspicuous oscillations. We have defined the ^He separation energy as 

53He = ^(Mg(Q)3He^3_i) - ^(Mg@3He^3) , (8) 

and have plotted this quantity as a function of N3 in Fig. [TOl S'sHe has a sawtooth structure 
similar to that of the atomic ionization energy vs. atomic number.— The more tightly bound 
^He atoms at A^3=8 and 20 allows one to identify shell closures, and correlate well with the 
local minima in Figs. |8] and [91 A similar result was obtained for Ca.-^ This is somewhat 
an unexpected result, as ^He droplets doped with Mg atoms look closer to axially than to 
spherically symmetric systems. We recall that for spin saturated systems, the first shell 
closures of the three-dimensional spherical harmonic oscillator appear for 2, 8, and 20 spin 
1/2 fermions. 



C. Absorption spectrum of Mg in^HejV4+^He7V3 clusters 

It is well-known that the shift and width of the electronic transitions of impurities in 
helium droplets are very sensitive to their environment and for this reason this technique 
is often employed to study their structure.-*^ We have calculated the dipole absorption 
spectrum of Mg as described in Ref. [l^ The line shape of the electronic transition is 
determined as 

i{uj) oc J dTz\^gs{n)\''s{oo + Vgs{n) - Kx(7^)) , (9) 

where {7^} refers to the positions of the atoms, and Vgs and Vex are, respectively, the ground 
and excited states potential energy surfaces. We recall that we have used the Mg-He X^T, 



interaction of Ref. 



14 for the 



24j for the ground state and the ^11 and potentials of Ref. 
excited states. To compute /(w) for a given value of u, one has to diagonalize a 3 x 3 matrix 
to determine the three components of the absorption line, each one arising from a different 



potentia' 
to Refs. 



energy surface, i.e., eigenvalue of the excited energy matrix. We refer the reader 
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and references therein for the details. 
The DMC calculation provides us with a set of walkers {TZj} representing the instanta- 
neous position of each atom in the cluster. These walkers have been used for determining 
the one-body densities presented before and can be also employed to obtain the absorption 
spectrum replacing / d7^|\E'gs(7?.)p by a sum over {TZj} in the above equation, so that a 
mixed estimator of I{uj) has been obtained. This is the same approximation employed for 
the calculation of the radial densities here studied finding no significant differences with 
their pure estimator values. 

As an illustrative example, the absorption spectrum of Mg is plotted in Fig. [TT]for three 
selected {N^, N4) combinations with + N4 = 40. The spectra are referred to that of Mg 
in the gas phase. As for other impurities, the long tail at high frequencies arises from the 
very repulsive contribution of the pair potential. The largest atomic shift corresponds 
to the ^He droplet and the smallest one to the ^He droplet, the one corresponding to the 
mixed droplet lies in between. 

The relative value of the shifts is easy to understand from the appearance of the helium 
densities shown in Fig. |2j Roughly speaking, the larger the density around the impurity, 
the larger the shift. For N3 + N4 = 40, Mg is fairly coated by helium and the value of 
the saturation density plays a significant role in the actual value of the shift. Since the 
saturation density is larger for ^He than for ^He this explains the relative position of the 
three absorption lines. 



IV. SUMMARY 

Using a DMC approach, we have found that a Mg atom in a '^He droplet experiences 
a transition from a surface to a bulk location for A4 > 26. For larger ^He droplets, the 
impurity resides in the bulk of the droplet. This conclusion has been drawn not only using 
the mixed estimator inherent to the importance sampling approach, but also using a pure 
estimator approach free from this bias. This finding agrees with the result obtained by Mella 
and coworkers, with which only minor quantitative differences are found. 

Due to the light mass of this alkaline earth atom and its weak interaction with helium, 
Mg impurities are found to be very delocalized inside droplets containing several thousand 
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^He atoms,— as those of interest in recent experiments.— In the case of ^He, we have 
found that Mg is always in the bulk of the ^He droplet, as determined by previous DFT 
calculations.^*^ To the best of our knowledge, the present calculations are the only DMC ones 
available for this system and for isotopically mixed helium droplets doped with Mg as well. 

For mixed droplets, we have found the well-known scenario^"^ that ^He and ^He atoms 
are distributed around the impurity into a onion-like shell structure, with ^He atoms coating 
the Mg impurity. This happens once the number of ^He atoms is large enough to fully cover 
the dopant. For small A^4 values, it may appear that some ^He is in contact with the Mg 
atom, which is again a consequence of the surface location of Mg in small ^He clusters. We 
want to point out, however, that distinguishing surface from volume regions in such small 
systems is largely arbitrary. 

The DMC walkers have been employed to calculate the dipole absorption line taking as 
a case of study three selected clusters with + = 40. We have found that the atomic 
shift is different enough for the studied configurations. This result only pertains to small 
mixed droplets. Indeed, when the number of ^He atoms in the droplet is large enough, the 
shift will be insensitive to the actual composition of the droplet, as the '^He component will 
be distributed into a shell distant from the impurity, whose absorption line will quickly tend 
to that of Mg in liquid '^ReM'^ 
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TABLE I: Parameters of the fit to the calculated Mg-He pair potentials. The values are such that 
Eq. ([2]) yields the pair potential in K when the distance r is expressed in A. 





X^T. (Ref. 24) 


(Ref. 14) 


Is (Ref. 14) 


A 


1.16902 X 10^ 


1.1308 X 10^ 


8.6195 X 10^ 


a 


3.0188 


1.6555 


0.00291 


/? 




0.17578 


0.14011 


D 


10.002 


9.3539 


21.075 


c& 


1.9454 X 10"^ 


2.0241 X 10^ 




Cs 


1.04717 X lO'^ 






Cio 




3.5046 X 10* 


2.3161 X 10^° 
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TABLE II: Ground state energies of Mg@ HeAr4 and Mg solvation energy (in K). 



N4 Ground state energy /XMg 








-15. yi 


(1) 


1 n on 
-lU.oU [2) 


on 


-5z.4o 


(2) 


-18. oz (7j 


21 


-55.99 


(2) 


-19.23 (5) 


22 


-59.35 


(2) 


-19.62 (b) 


26 


-Oz.o4 


(2) 


-19.00 [1) 


24 


-00.44 


(2) 


on 00 /'n\ 
-20.33 (9j 


zo 


-7U.U7 


(2) 


-20.57 (8j 


26 


-73.75 


(2) 


-20.99 (9) 


27 


-77.46 


(2) 


-21.46 (9) 


28 


-81.13 


(2) 


-21.89 (9) 


29 


-84.87 


(2) 


-22.17 (7) 


30 


-88.74 


(2) 


-22.63 (10) 


40 


-128.96 


(3) 


-25.83 (12) 
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TABLE III: Ground state energies of Mg@^Hejv4+^Hejv3 and Mg solvation energy (in K). 



AT AT 


Ground state energy 


MMg 


z z 


-b.olz [0) 


-o.zOl (oj 


8 8 


-28.77 (2) 


-16.89 (2) 


20 8 


-69.73 (4) 


-23.03 (4) 


8 20 


-43.13 (2) 


-22.62 (3) 


20 20 


-88.51 (5) 


-26.25 (6) 
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TABLE IV: Ground state energies of MgO^HeATg (in K). 





Ground state energy 


Ns 


Ground state energy 


2 


-1.887 (2) 


13 


-9.31 (1) 


3 


-2.544 (3) 


14 


-9.67 (1) 


4 


-3.312 (3) 


15 


-10.64 (1) 


5 


-4.206 (4) 


16 


-11.61 (1) 


6 


-5.124 (4) 


17 


-12.76 (1) 


7 


-6.141 (5) 


18 


-13.89 (1) 


8 


-7.266 (5) 


19 


-15.19 (1) 


9 


-7.666 (6) 


20 


-16.45 (1) 


10 


-8.055 (7) 


21 


-16.820 (1) 


11 


-8.508 (7) 


40 


-33.57 (3) 


12 


-8.885 (8) 
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TABLE V: Comparison of several DMC results for Mg@^HeAr4 clusters. 





Cluster energy (K) 


Mg solvation eneri 




AT 
iV4 


inis work 


Ker. JA Ker. ^ 


Inis work 


Kei. .14 


of 09 


o 
8 


-15.91 (1) 


-15.91 (1) 


-10.80 (2) 




-10.76 (i) 


20 


-52.46 (2) 


-51.64 (1) 


-18.52 (7) 


-18.49 (1) 




25 


-70.07 (2) 


-69.04 (1) -70.15 (6) 


-20.57 (8) 


-25.03 (4) 


-20.93 (9) 


30 


-88.74 (2) 


-87.41 (1) 


-22.63 (10) 


-26.89 (4) 




40 


-128.96 (3) 


-126.51 (3) 


-25.83 (12) 


-25.39 (4) 
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FIG. 1: (Color online) Pair potentials used in this work. The symbols correspond to the results 
of Hinde2^ for the ground state (left panel) and Mella et al}'^ for the excited states (right panel). 
The lines correspond to the parameterizations given in Table H The Aziz He-He potential^^ has 
also been plotted in the left panel (solid line). 
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FIG. 3: (Color online) Atomic density of the mixed drop with A'4 = 8, = 20. The contour plots 
of ^He and ^He densities are separately displayed at the top of the figure. The darker the region, 
the higher the density. The three-dimensional plot of the same densities is displayed at the bottom 
of the figure, with the same plot box as in the top contours. 
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FIG. 5: Root mean square radii of the He clusters and the Mg impurity, measured with respect 
to the center-of-mass of the He cluster in both cases. The dotted lines are to guide the eye. 
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FIG. 6: Grey scale two-dimensional plot of the particle density (right panel) for different 
Mg@^Hejv4 clusters; the lighter the region, the higher the helium density. The He densities are 
displayed in boxes of 30 A x 30 A. The left panels represent the probability distribution P(r) of 
finding the Mg atom at a distance r from the He center-of-mass [/q°° drP{r) = 1]. 
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FIG. 7: (Color online) Evolution of the rms radii with the length of blocks of the pure estimator. 
The dashed lines are to guide the eye. 
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FIG. 8: (Color online) Ground state total energies (in K) of MgQ^HeAr^+'^HeAra droplets. The 
dotted lines are to guide the eye. 
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FIG. 9: (Color online) Solvation energy of Mg (in K) in MgQ^HeAr^+^HeATg droplets. The dotted 
lines are to guide the eye. 
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FIG. 10: Separation energy of a ^He atom (in K) in MgQ^HeATg droplets. The dotted line is to 
guide the eye. 
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FIG. 11: (Color online) Dipole absorption lines of Mg in three selected N3 + = 40 droplets. 
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